Subtraction of temperature induced phase noise in the LISA frequency band 
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Temperature fluctuations are expected to be one of the limiting factors for gravitational wave 
detectors in the very low frequency range. Here we report the characterisation of this noise source 
in the LISA Pathfinder optical bench and propose a method to remove its contribution from the 
data. Our results show that temperature fluctuations are indeed limiting our measurement below 
one millihertz, and that their subtraction leads to a factor 5.6 (15 dB) reduction in the noise level 
at the lower end of the LISA measurement band (10 -4 Hz), which increases to 20.2 (26 dB) at even 
lower frequencies, i.e., 1.5 x 10 -5 Hz. The method presented here can be applied to the subtraction 
of other noise sources in gravitational wave detectors in the general situation where multiple sensors 
are used to characterise the noise source. 



I. INTRODUCTION 

Temperature fluctuations are expected to be one of 
the limiting noise contributions to gravitational wave in- 
terferometers at low frequencies, in the millihertz band. 
Currently, ground-based gravitational wave detectors are 
not limited by this low frequency contribution because 
they are already dominated by stronger noise contribu- 
tions at low frequencies [1 j, primarily due to seismic noise 
which usually encompasses a variety of effects from hu- 
man activity, to environmental microgravity effects. This 
is one of the main driving forces behind the proposal 
of LISA [2 j, a space-borne gravitational wave detector 
which aims at observing the gravitational wave sky in 
the low frequency region of the spectrum, meaning fre- 
quencies down to f ^ O.lmHz. The LISA concept is a 
constellation of three spacecrafts, each one located at a 
vertex of a 5 million kilometers triangle. The constella- 
tion follows the Earth in a heliocentric orbit, 20° behind 
the Earth, with a 1° inclination with respect the eclip- 
tic. Each spacecraft hosts two test masses in nominally 
geodesic motion. Inter-spacecraft laser interferometry is 
used to measure differential displacements between these 
test masses, which contains the gravitational wave sig- 
nal. Although ESA and NASA ended their collaboration 
pursuing the implementation of this mission, ESA is cur- 
rently considering a mission which will gather much of 
the existing LISA heritage, the eLISA mission [3], known 
within ESA as NGO (New Gravitational wave Observa- 
tory). 

For any low frequency space-borne gravitational wave 
detector, the main contributions limiting the instrument 
performance at low frequencies are expected to come 
from spurious accelerations. Thermal induced effects are 
expected to play an important role in these contribu- 
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tions. Although thermal isolation in the spacecraft is ex- 
pected to account for a ~99% reduction of the remaining 
solar temperature modulations [4 , the instrument sen- 
sitivity will be directly coupled to thermal effects pro- 
ducing forces on the test masses as is the case of the ra- 
diometer effect, radiation pressure, outgassing [5 or, as 
lately outlined, brownian gas motion [6]. Also, thermal 
related distortions of the optical path or the influence of 
temperature on other interferometer components, such 
as reference cavities, modulators or photodiodes could 
have an impact on the instrument performance. The 
ability to prevent or, if possible, remove any effect at 
the low frequency end of the LISA sensitivity curve is 
highly desirable since the low frequencies could contain 
valuable information about black-hole binaries resulting 
from mergers of pre-galactic structures and galaxies [7]. 

LISA Pathfinder [HJ [9] is an ESA mission, with some 
NASA contributions, that will test key technologies re- 
quired for LISA. In particular, it will explore and char- 
acterise these thermal related effects around 1 mHz. 

Ground experiments to test LISA technologies and 
LISA Pathfinder flight hardware are the first ones fac- 
ing this problem and as such, they are the natural play- 
ground to investigate this noise source. Below millihertz 
frequencies, daily temperature modulations enter into the 
measurement band and are poorly screened since the re- 
quired isolator for that would need a high amount mass, 
assuming a passive insulator [10]. A different approach 
is to actively control the thermal fluctuations at low fre- 
quencies [11 , or in other words, to correct the slowly 
varying thermal drift. However, this shifts the problem 
to the design of a highly stable control loop able to mea- 
sure and remove the low frequency fluctuation with high 
precision. 

Here we present an alternative to isolation which is to 
measure the temperature and subtract its contribution 
from the main interferometric data stream. Our method 
computes the temperature contribution to interferometer 
data by first determining the transfer function between 
both measurements. This is then translated into a digi- 
tal filter which is used to compute the noise contribution 
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FIG. 1. Schematic representation (not to scale) of the loca- 
tion of the temperature sensors in the LISA Pathfinder optical 
bench. One sensor was attached to each mounting structure 
supporting the mirrors, a third sensor was attached to the 
side slab flange — a metallic structure from which the optical 
bench will be supported inside the satellite thermal shield — 
and a fourth (not shown) measured the temperature in the 
laboratory. 



quired drag-free technology used in the final mission is 
part of our setup. Although the test-bed has shown the 
required performance in the LISA Pathfinder frequency 
band, lmHz < f < 1 Hz [I3j [14], our main aim was 
to determine to what extent our measurements were lim- 
ited by environmental temperature fluctuations and, also, 
to evaluate this noise contribution in the LISA band. 
In order to do that, we measured the temperature in 
four different locations of our setup, as shown in Fig. [T] 
We initially tested the method using a long data set 
(2 x 10 5 s), where the interferometer was running without 
active noise suppression control loops. With the temper- 
ature coupling characterised and the method verified, we 
then proceeded to subtract the temperature contribution 
in a shorter, low noise data segment to characterise the 
impact of this contribution when the interferometer is 
performing according to the mission requirements, reach- 
ing 1 pm/\/Hz at 0.1 Hz. 



III. EXPERIMENT CHARACTERISATION 
A. Signal processing definitions 



that is finally subtracted from the main interferometric 
measurement. A crucial point in this scheme is that, in 
evaluating the phase of the transfer function, we take 
into account the group delay between the time of the 
temperature measurement and the actual effect in the 
interferometer so that the subtraction is performed co- 
herently. 

The method is general and can be applied to the sub- 
traction of any noise contribution with a delayed impact 
on the measurement. To do this, the variable driving 
the noise contribution (temperature in our case) must be 
measured. As we show below, the error assigned to the 
subtraction process will be computed from the coherence 
between this magnitude and the data once this is sub- 
tracted. 

This paper is structured as follows. In Sec. ITT] we de- 



scribe our experiment setup, in Sec. [HI] we introduce the 
notation and the basic definitions that we then use to 
characterise our system. In Sec. IV we propose two 
analysis schemes to subtract the temperature noise con- 
tributions from the interferometer, which we then apply 
in Sec. [V] We discuss the results and conclude in Sec. lVIl 



II. SETUP DESCRIPTION 

The measurements used in our analysis were taken 
in the LISA Pathfinder engineering model optical 
bench [12] , a close replica of the final bench to be flown in 
the LISA Pathfinder satellite. In our setup, test masses 
are substituted by fixed mirrors so only the interfero- 
metric metrology subsystem is tested; none of the re- 



We introduce here the basic notation needed to de- 
scribe our experiment. In general, we could describe our 
experimental data as a set of q measured temperature 
time series, T^(t), i — 1 . . . which pass through q sys- 
tems and combine to produce a single measured inter- 
ferometer output, $(£). If we further assume that our q 
systems are linear systems, this can be directly translated 
into the frequency domain as a system of equations 



$(u)) = ^H i (u))T i (u)) + n(u)) 



(1) 



where <l>(co>) and T{uS) are the Fourier transforms of the 
input and the output respectively, and H(uo) is the sys- 
tem transfer function. The transfer function is usually 
estimated from real data by 



(2) 



where St^^) and 5^7; (<^0 are the cross-power spectral 
density and the power spectral density, which are defined 
as the Fourier transform of the correlation functions. A 
second quantity relevant in our case is the coherence func- 
tion, 
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(3) 



which quantifies the amount of correlation between two 
data sets in each frequency bin. The coherence function 
turns out to have a particular relevance in our analysis 
since it allows us to quantify the error on the estimation 
of the transfer function. To bound the uncertainty of our 
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FIG. 2. Setup characterisation. We show the matrix of coherence functions between the different measurements. In the upper 
raw the coherence between interferometer phase and temperature sensors, the rest of figures provide the coherence functions 
between the different temperature sensors. Units for the x's values are in hertz. 



transfer function estimates, we compute the error on the 
transfer function as [15] 

a [\H Ti *\] * a [0(H Ti *)] ~ & " 7 f'*f^ , (4) 

where rid is the number of averages used to compute the 
transfer function estimates. The error on the transfer 
function already contains the information about the cor- 
relation between the temperature and the interferometer 
in the form of the coherence function 7^ T . (/). By doing 
this, we ensure that our error estimate takes into account 
the correlation between temperature and interferometer 
phase fluctuations. Notice that the same applies to the 
correlation between temperatures, where we will be using 
the coherence 7^. T .( CJ ) t° evaluate the interdependence 
between temperature variations at each sensor. 



B. Data preprocessing and characterisation 

Before analysing, the data needed to be resampled onto 
a common time grid. Both the interferometer data (with 
a sampling frequency of f s = 32.47 Hz), and the tem- 
perature data (with f s = 1.3 Hz) were down-sampled to 



f s = 1 Hz. The temperature sensor data were interpo- 
lated to the new time grid with the new sampling fre- 
quency, from 1.3 to 1 Hz. In the interferometer case, the 
data are first down-sampled by a factor of 10 in order to 
ease data handling, and then resampled to the common 
grid. 

To determine the interdependence between tempera- 
ture and interferometer readout, we compute the coher- 
ence functions between them . As shown in the up- 
per raw of Fig. [2j the three sensors attached to the bench 
show a coherence above 80% below / ~ 0.3 mHz. As ex- 
pected, the sensor measuring temperature variations in 
the lab has the lowest coherence with the interferometer, 
although still a coherence of 80% can be assigned between 
interferometer and temperature variations in the lab en- 
vironment. While at higher frequencies environmental 
temperature fluctuations can be screened by the vacuum 
tank; in the low frequency region these fluctuations pass 
directly through the setup (although delayed, as we show 
below). 

We proceed by determining the interdependence be- 
tween temperature measured at different locations. In 
this case, the coherence function is JtiTj (f)U^j- Results 
in Fig. [2] show that the four temperatures are strongly 
correlated (7^.7^ (/) — 1)? m the vei T i° w frequency fre- 
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FIG. 3. Temperature to interferometer transfer functions for 
the four sensors in the setup. Errors and fitted model are only 
shown for the transfer function that shows the highest con- 
tribution, the one coming from the sensor in the first mirror 
(Ml). 



quency band (/ < 0.5mHz), and the correlation drops 
off rapidly above these frequencies. Hence, at the low 
frequencies of interest, the four measurements contain 
mostly the same information regarding coupling between 
interferometer phase and temperature. 

Even though the four measurements provide similar 
temperature information at low frequencies, they do not 
contribute with the same strength to the interferometric 
measurement. To show that, we use the transfer func- 
tion between the temperature data and interferometer 
data -Eq. In Fig. [3] we show these transfer func- 

tions for each sensor; the sensor attached to the first test 
mass (Ml) shows a higher contribution with respect to 
the other — numerical values for the transfer function are 
shown in Tab.|IJ Also, from a linear fit to the phase of this 
transfer function in the low frequency band we obtain the 



group delay of the filter. For our setup, the response of 
the interferometer to temperature fluctuations below the 
millihertz frequency band is delayed by 1430 ± 70 s with 
respect to the actual time of temperature measurement. 

We note that, in our analysis, and as shown in Tab. [l| 
this correlation is valid for f < 0.5 mHz. Therefore, the 
frequency region where we would expect a noise reduction 
after subtraction of this contribution would be below the 
LISA Pathfinder measuring band, but directly affecting 
a LISA-like measurement. 



IV. SPECTRAL ANALYSIS 

Once we have characterised our system, our next step 
is to disentangle the temperature contribution from the 
interferometer measurement. We will proceed in two dif- 
ferent ways: firstly, by solving the linear system and ob- 
taining a set of optimal transfer functions and secondly, 
by applying a conditioned scheme, which proceeds by 
solving the linear system for each individual contribu- 
tion. We describe both in the following section. 

A. Optimal spectral analysis 

In order to obtain the system of linear equations that 
will describe our problem we start with Eq. 0, 

q 

= ^Hi{u)Ti[u)) ± N(l>), (5) 

where, if we multiply &(uj) by its complex conjugate 
<£(cj)*, and taking expected values we obtain 

q q q q 

S nn = Hi S^n — ^ Hj ftr-s+y^ Hj Hi St^Ti • 

i=l j = l i=j j=l 

(6) 

The optimal solution for this system of equations which 
minimise S nn over all possible choices of Hi is found by 



TABLE I. Numerical values for the transfer function corre- 
sponding to the sensor in the first mirror mount (Ml) and 
associated errors. Last column shows the coherence function. 



Frequency [Hz] 


\H vT \ [rad/K] 




1.6 x 10" 5 


4.5 ±0.5 


0.97 


5.1 x 10" 5 


5.1 ±0.6 


0.96 


1.1 x 10" 4 


4.5 ±0.6 


0.91 


2.6 x 10~ 4 


2.7 ±0.5 


0.85 


5.3 x 10" 4 


1.6 ±2.8 


0.12 


1.0 x 10" 3 


1.4 ±2.6 


0.10 


2.5 x 10" 3 


0.5 ±6.5 


0.01 
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FIG. 4. Comparison between optimal and conditioned analy- 
sis. We show in the same figure the LISA Pathfinder require- 
ment as a reference. 
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dHj 

which leads to f!5l 



oh* 



o, 



(7) 



5>s 



j = 1,2, ...q, 



(8) 



which can be solved for i^, since SW* and S^Tj are 
known. In our particular case this turns into a system 
of q = 4 equations, which leads to analytical solutions 
for Hi, i = Deriving the expressions for 

which we do not reproduce here, requires us to compute 
24 terms for each transfer function, each term containing 
a product of 4 spectral densities. In the following section 
we introduce a sequential method which allows us to han- 
dle simpler equations, leading to equivalent results. 



B. Conditioned spectral analysis 

A second possible approach is to subtract (in a linear 
least square sense) one by one those contributions per- 
turbing our measurement. This scheme leads to simpler 
expressions, which will allow the derivation of the digi- 
tal filters required to clean the data in the time domain. 
This approach, known as conditioned spectral analysis, 
derives the dependence between one of the inputs and 
the output when the other inputs are turned off, implic- 
itly requiring here that the correlations between inputs 
are turned off as well. In order to introduce this analysis 
scheme, we refer back to Eq. which we can write now 



= J2L Ti ^oj)T Hi _ lV Xoo) +N(u>), 



(9) 



i=i 



where T i .^ i _ 1 y(oj) represents the Fourier transform of 
variable Ti(t) when the linear effects of T\ to T^_i have 
been sequentially removed from Ti by optimum linear 
least-squares techniques. Notice that each of this or- 
dered conditioned records will be mutually uncorrelated, 
a property which is not generally satisfied by the origi- 
nal records. The frequency response function L$Ti is the 
optimum linear system to predict <fr(t) from Ti(t). Anal- 
ogously as we did in Eq.([6|, we can derive the optimum 
operator, defined as the one that minimises S nn for any 
possible combination of L^tv This leads to [15] 



^Ti<S>-T {i _ iy . 



s> 



T i T i -T (i _ 1)l 



(10) 



where the spectra and cross-spectra appearing in the pre- 
vious expressions are computed on the conditioned vari- 
ables. Notice that for the case of a system where q = 1, 
Eq. (10) and Eq. ([8| will naturally reduce to the same 
expression 



= Hi = 



Sri 



(11) 



Following Eq. ([3| we can define the coherence function 
for conditioned variables as 



(i-Di s- 



T l T z -T (i _ 1)] O^.T (i _ 1)] 



(12) 



i.e., the coherence between <1> and T^, once all contribu- 
tion from the previous i — 1 sensors, T^_iy^ are sub- 
tracted. Together with Eq. [lOj it can be shown that the 
conditioned spectra can be written down in terms of par- 
tial coherence functions as follows [15] 



SW-Tii = SW-T^x), (l 



(i-iy. 



(13) 



where i = 1,2, ...g. This quantity corresponds to the 
phase noise spectrum, S$$(/)j when contributions from 
temperature sensors, Ti are recursively removed by linear 
least squares. 



We will use Eq. ( 13 ) in the following to obtain predic- 



tions of the expected noise reduction when we subtract 
a given temperature contribution from the main interfer- 
ometer measurement. This will give us a useful diagnos- 
tic tool to evaluate the contribution coming from each 
location independently. Also, it is worth stressing that, 
in comparison with the optimal method, the sequential 
approach in the conditioned scheme allows us to deal with 
expressions with a maximum of q terms, instead of the 
q 2 linear system required in the optimal approach, being 
q the number of noise contributions in our problem. 
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FIG. 5. Temperature noise subtraction for a two days measurement without control loops suppressing the interferometric noise 
contributions. Left: Time domain data series for the two most relevant temperature sensors and the interferometer read out. 
We also show interferometer data after subtracting the information contained in the two previous temperatures sensors. Right: 
Comparison of the frequency domain analysis — previously compared with the optimal prediction in Fig. [4] — with the spectra 
of the interferometer time series after the subtraction of the temperature contribution coming from the two most relevant 
temperature sensors. For comparison, we also show the requirement on interferometric noise in the LISA Pathfinder mission. 



V. TEMPERATURE NOISE SUBTRACTION 

A. Frequency domain analysis 

Before removing the contribution arising from temper- 
ature, we proceed to estimate the expected reduction. To 
do this, we compute the conditioned spectra. Results are 
shown in Fig. [4] where the curve labelled as 'prediction' 
shows the conditioned spectrum in Eq. (13) for the two 
sensors with the strongest contributions — according to 
Fig. [3] We also performed the analysis for a third sensor 
but the coherence with the interferometer was already 
too small to show any improvement. In Fig. [4] we also 
compare the conditioned spectra for the two first subtrac- 
tions with the spectrum obtained following the optimal 



subtraction scheme, presented in Sec. IV A The compar- 
ison confirms that the second subtraction in the condi- 
tioned scheme achieves the optimal level and, therefore, 
subsequent subtractions will not improve substantially. 
As expected, a coherent subtraction of the temperature 
contribution would reduce the noise level at frequencies 
f < 0.4 mHz. According to this first analysis, the sub- 
traction of the contribution contained in these two sen- 
sors would reduce the instrument noise floor by a factor 
5.6 (15 dB) at the very low end of the LISA measure- 
ment band, f = 0.1 mHz. The factor increases up to 20.2 
(26 dB) at the lowest frequency bin, f = 15 /iHz. 



B. Time domain analysis 

We want to be able to subtract noise contributions in 
time domain, which can be useful to avoid effects purely 
related to the Fourier domain, e.g. correlation between 



frequency bins. Time domain analysis is of relevance as 
well in order to obtain temperature noise cleaned time- 
series which can be used for subsequent analysis. For 
these reasons we have previously resampled both data 
streams with the same sampling frequency and on to a 
synchronous time grid. As previously shown, the relation 
between temperature and phase read out is better de- 
scribed by their frequency domain description. We then 
translate this transfer function into a digital filter, i.e., a 
recursive relation that allows us to include a certain dy- 
namical response and, in particular, a delayed action of 
the temperature upon the interferometer. The tool used 
here is the vector fitting algorithm [16] which allows us 
to fit the measured transfer function in terms of N poles, 
Pk : and residues , Zk, 



N 



k=l 



■Pk z~ 



(14) 



Since our previous analysis showed that the tempera- 
ture contribution to the interferometer is relevant in a fre- 
quency range up to ~ 1 mHz, we fit the transfer function 
up to this frequency. Figure [3] shows the fit result with a 
3rd order model. Once the digital filter is obtained, we 
filter the temperature measurement with /i<£T;(£ -1 ) to 
obtain the temperature contribution to the interferome- 
ter, which we can then readily subtract from the original 
measurement. This procedure is performed initially con- 
sidering the temperature reading of the sensor attached 
to the mounting of the first mirror (Tmi)- Then, the in- 
formation contained in the flange sensor (Tfl) is removed 
from the residual of the first subtraction. 

As shown in Fig. [5j the subtracted curve is in agree- 
ment with the one previously obtained in the frequency 
domain analysis. We show, for comparison, the LISA 
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FIG. 6. Temperature noise substraction for a data segment 
where control loops were active, showing compliance with the 
Pathfinder interferometric noise requirement. 



Pathfinder interferometric noise goal on the same plot. 
Since the control loops were not suppressing the interfer- 
ometric noise contributions, the measurement is above 
the goal. To further investigate the temperature noise 
contribution in our setup we used the method above to 
subtract the temperature noise contribution in a mea- 
surement where the instrument was operating in closed- 
loop (all noise suppression on). The transfer function is 
the same as the one previously evaluated since the setup 
remains unchanged. As shown in Fig. [6j the effect of the 
subtraction reduces the noise level at lower frequencies, 
f < 0.2 mHz, when compared with the control- free case 
although the effect is stronger at f = 0.1 mHz, reducing 
the noise floor a factor 10 (20 dB) after the subtraction. 
We notice as well that, according to our current analy- 
sis, we can not attribute the complete noise contribution 
in the low end of the LISA Pathfinder measuring band, 
f ~ 1 mHz, to temperature driven phase fluctuations. 



VI. DISCUSSION AND CONCLUSIONS 

The analysis reported here characterises, for the first 
time, the low frequency temperature coupling to an inter- 
ferometric measurement in a realistic gravitational wave 
instrumental setup. We also introduce a methodology to 
clean the main measurement of this noise contribution by 
removing the information contained in multiple sensors 
monitoring the experiment. Once applied, we obtain a 
noise-free time series whose spectral density confirms the 
prediction obtained by a frequency domain noise projec- 
tion. We have shown as well that the proposed method, 
based on conditioned data streams, agrees with the pre- 
diction obtained with an optimal approach. The advan- 
tage of the former method being a simplification in the re- 
quired analysis as well as a straightforward characteriza- 
tion of the contribution of each temperature sensor inde- 
pendently. Results on the LISA Pathfinder optical bench 
show that, when considering the temperature readings 
in our current setup, the instrument performance is not 
affected in the Pathfinder measurement band after the 
temperature noise subtraction, although it gives a signif- 
icant noise reduction at lower frequencies, and therefore 
this noise source has a direct impact for LISA-like experi- 
ments. This is of relevance for any ground-based gravita- 
tional wave experiments aiming at low frequencies, and in 
particular for those testing technologies for space-based 
gravitational wave observatories. It is worth stressing 
that, after temperature noise subtraction, our experi- 
ment shows a performance of 100nm/A/Hz at f = 30 
/iHz. This precision at very low frequencies in an optical 
metrology measurement can only be compared with re- 
sults from a recent test campaign of the LISA Pathfinder 
spacecraft in an artificial space environment [17], which 
further confirms that the method described here can al- 
low on-ground lab experiments to compute a performance 
curve without the limitation of the unavoidable temper- 
ature drifts, achieving results which are close to space 
conditions. The method is extendable as well to any 
frequency dependent contribution that is independently 
measurable in gravitational wave experiments. 
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